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Abstract 

We present Hanbury-Brown Twiss measurements from the PHENIX experiment at RHIC for 
final results for charged kaon pairs from y^AW = 200 GeV Au+Au collisions and preliminary 
results for charged pion pairs from ^fs = 200 GeV p+p collisions. We find that for kaon pairs 
from Au+Au, each traditional 3D Gaussian radius shows approximately the same linear increase 
as a function of A^ r 3 t . An imaging analysis reveals a significant non-Gaussian tail for r > 10 fm. 
The presence of a tail for kaon pairs demonstrates that similar non-Gaussian tails observed in 
earlier pion measurements cannot be fully explained by decays of long-lived resonances. The 
preliminary analysis of pions from s/s - 200 GeV p+p minimum biased collisions show cor- 
relations which are well suited to traditional 3D HBT radii extraction via the Bowler- Sinyukov 
method, and we present R 0VLt , R S id Q , and R\ ong as a function of mean transverse pair mass. 



1. Introduction 

Measurements of the Hanbury-Brown Twiss effect for hadrons in relativistic heavy-ion col- 
lisions provide information on the source size and hadron emission duration, allowing a map- 
ping of the freeze-out hyper- surface. Advancement of both HBT experimental analysis methods 
and hydrodynamic modeling must continue in order to extract precision information, and recent 
progress in resolving the historic discrepancies between HBT measurements and hydrodynamic 
model results Q emphasizes this point. The observations of extended, non-Gaussian, source 
size from two-pion correlations [2] make the measurement of two-kaon correlations important 
for understanding the contribution from decays of long-lived resonances. 

HBT measurements show azimuthal sensitivity relative to the reaction plane of non-central 
collisions. Future femtoscopic measurements may similarly use a jet axis to define the event- 
by-event geometry to gain insight into the modification of jets by the Quark-Gluon Plasma, 
and conversely, the impact of the jet on the medium. Such measurements will require a better 
understanding of HBT in p+p collisions. 

The main source of information for HBT analyses is the two particle correlation function 
C2(q) where, for identical bosons, Bose-Einstein enhancement is manifest as an increased yield 
at low relative pair momentum q. Several forms for C2(q) are shown in Eq[T] The first form is 

A(q)/B(q) = C 2 (q) * 1 + J dr(\%(r)\ 2 - 1)5 (r) * 1 - X + AF c (q)[l + G] (1) 

the experimental correlation function where A(q) is the measured two particle pair momentum 
difference distribution and B(q) is the normalized background distribution. In the second form, 
the Koonin-Pratt equation, r is pair separation, O q (r) is the relative wave function and S(r) is 
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the emission source. In the final form, Bowler-Sinyukov, F c (q) is the Coulomb correlation and 
G = exp(-^ ide ^ ide - R 2 0Ut q 2 0Ut - ^ 2 ong ^ 2 ong ). We discuss this in more detail below. 

2. Methodology 

Charged kaons and pions are reconstructed using the PHENIX Drift Chambers (DC), Pad 
Chambers (PC), and lead scintillator Electromagnetic Calorimeters (EMCal). The DC and PC 
track particles through a magnetic field to obtain momentum, and the EMCal provides time of 
flight relative to the Beam-Beam Counters (BBC). Matching cuts are applied between detector 
hits and track projections at the EMCal and DC farthest from the vertex. This information is used 
to calculate a mass-squared for particle identification. The combination of BBC and Zero-Degree 
Calorimeters provide vertex and centrality information. 

The experimental correlation function in Eq.[T]is obtained from the ratio of the relative mo- 
mentum distribution of all pairs within same events, A(q), to the normalized relative momentum 
distribution of pairs from different events within the same event class, B(q). Event class refers 
to events with similar centrality, collision vertex et cetera. Pair selection cuts are used to remove 
track splitting and merging. For the analysis of kaons from Au+Au, Monte-Carlo detector sim- 
ulations are used to correct for reconstruction inefficiencies for particle pairs with close paths 
through the DC and EMCal. These inefficient regions are removed by pair cuts for preliminary 
analysis of pions from p+p. 

Information about the emission source is obtained from C2(q) using two methods. The tra- 
ditional approach makes assumptions about the emission source to derive a form for C2 with pa- 
rameters which we fit using a log-likelyhood method to best match the measured A(q) and B(q). 
The Bowler-Sinyukov formulation, which assumes a Gaussian emission source at the core and a 
long range halo, has become the standard formulation of this approach. Using Bowler-Sinyukov 
requires the choice a coordinate system and reference frame. We choose the out-side-long coor- 
dinate system, where long is the beam axis direction, out is the direction of the pair transverse 
momentum and side is orthogonal to both. The Longitudinally Co-moving System (LCMS), 
where total longitudinal pair momentum is zero, is chosen as the reference frame. The second 
approach is to decompose the measured C2 - 1 using suitable basis functions, and then directly 
solving the integral equation for S (r) (3). This process, referred to as imaging, requires no strong 
assumptions about the form of S(r), and no explicit Coulomb correction is needed since the full 
wave function is used. 

3. Charged kaons from Au+Au collisions 

We refer the reader to |4) for additional analysis details and a full presentation of the data; 
however, the centrality dependence of the Gaussian source parameters is shown in Fig. [T] Also 
shown are the pion data from the same system rescaled to the transverse mass of the kaon data. 
The rescaled pion data overlap with kaon data, and both show a linear trend. A simultaneous fit 
of all kaon radii to a line with zero intercept produces a slope of 0.51 ± 0.01 and a^ 2 /DoF 
of 14.6/8. The mr dependence of the Gaussian source parameters shows poor agreement with 
traditional hydrodynamic models which match flow and spectra measurements |4 |. A recent ID 
hydrodynamic calculation [T], which improves on several aspects of traditional hydrodynamic 
models and includes novel universal initial flow, reproduces the kaon radii well even though it 
was not tuned for this data set. The model must now be extended to higher dimensions to test 
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Figure 1: Centrality dependence of Bowler-Sinyukov source parameters for charged kaon pairs compared to pion data 
1 5 1 rescaled to the same transverse mass. The dashed line is the result of a simultaneous fit of all kaon radii to y = mx. 




Figure 2: Slices of the 3D correlation function for n + pairs. Solid curves show results for several fit forms of C2. 



if it can simultaneously describe the flow data. A ID imaging analysis was also performed l5l . 
The reconstructed source shows Gaussian emission at low r, but a significant non-Gaussian tail 
is observed for r > 10 fm. 

4. Charged pions from p+p collisions 

Roughly 2.5 million like sign pion pairs from minimum biased p+p collisions from the 2004- 
2005 RHIC run are analyzed. Despite concerns regarding the smoothness approximation [6] 
and energy-momentum correlations for small and low multiplicity systems 0, many ID and 
2D correlation slices from p+p collisions were examined and verified to be well behaved and 
sufficiently described by the Bowler-Sinyukov form for our statistics and acceptance. Example 
ID slices of the full 3D correlation are shown in Fig. [2] The Bowler-Sinyukov method is used 
to extract the emission source parameters show in Fig. [3] We note a correction from the oral 
presentation, which swapped R out and R S ^ Q labels. 
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Figure 3: The transverse mass dependence of the Bowler-Sinyukov source parameters for pion pairs from = 200 GeV 
minimum biased collisions. We note a correction from the oral presentation, which swapped R out and R^e labels. 



5. Discussion and Outlook 

We find that for kaon pairs from Au+Au, R ouU R S id&, and R\ ong show approximately the same 
linear dependence of R « (0.51 ± 0.001)A^ r 3 t . An imaging analysis shows that the bulk emission 
at low radius is consistent with a Gaussian source but a significant non-Gaussian tail is revealed 
for r > 10 fm. The presence of this tail for kaon pairs indicates that similar non-Gaussian 
tails observed in earlier pion measurements have a component other than decays of long-lived 
resonances, greatly increasing the likelihood of a truly extended emission source. 

The preliminary analysis of pions from yfs = 200 GeV p+p minimum biased collisions 
show correlations which are surprisingly well suited to traditional 3D HBT radii extraction via 
the Bowler-Sinyukov method for our statistics and acceptance. An imaging analysis to provide 
further insight will be performed in the near future. 

The next era of HBT should be an exciting time of learning from hydrodynamic models that 
simultaneously describe spectra, flow, and HBT as well as experimentally exploring correlations 
relative to the jet axis in p+p and, eventually, heavy ion collisions. 
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